Baryon Number Fluctuation and the Quark-Gluon Plasma 
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We show that lob or lo b , the squared baryon or antibaryon number fluctuation per baryon or 
antibaryon, is a possible signature for the quark-gluon plasma that is expected to be created in 
relativistic heavy ion collisions, as it is a factor of three smaller than in an equilibrated hadronic 
matter due to the fractional baryon number of quarks. Using kinetic equations with exact baryon 
number conservation, we find that their values in an equilibrated matter are half of those expected 
from a Poisson distribution. Effects due to finite acceptance and non-zero net baryon number are 
also studied. 



I. INTRODUCTION 

A new state of matter, the quark-gluon plasma, is 
expected to be formed in heavy ion collisions at ultra- 
relativistic energies, such as at the Relativistic Heavy 
Ion Collider (RHIC) that has just begun its operation at 
the Brookhaven National Laboratory. Many observables 
have been proposed as possible signatures for the quark- 
gluon plasma phase during the collisions pL such as 
strangeness enhancement J/ip suppression ||, mod- 
ification of high px particle spectrum [Q, and Mj- scal- 
ing H and double phi peaks || in the dilepton spec- 
trum. Recently, event-by-event fluctuations of various 
particles have also attracted much attention 0. Since 
the baryon and charge numbers of quarks are fractional, 
the net baryon number and net charge || fluctuations 
per particle are smaller in the quark-gluon plasma than 
in the hadronic matter and have also been suggested as 
possible signatures for the quark-gluon plasma. 

Because of conservation of total net baryon number, 
the net baryon fluctuation is finite only if one considers a 
part of the full phase space, e.g., in a finite rapidity range. 
Since the quark-gluon plasma is eventually converted to 
the hadronic matter, elastic scatterings of baryons and 
antibaryons would modify their phase space correlations 
and thus increase the net baryon fluctuation towards that 



expected from an equilibrated hadronic matter |10|. On 
the other hand, elastic scatterings of baryons and an- 
tibaryons do not change their numbers in the full phase 
space, so fluctuations of baryon and antibaryon numbers 
are alternative signatures for the quark-gluon plasma in 
relativistic heavy ion collisions. 

In this study, we shall use a master equation with ex- 
act baryon number conservation to derive the baryon and 
antibaryon number fluctuations in an equilibrated mat- 
ter, and to study the effects due to finite acceptance and 
non-zero net baryon number. 



II. KINETIC MODEL 

Due to the different fundamental units of baryon num- 
bers in quark-gluon plasma and in hadronic matter, fluc- 
tuations of baryon and antibaryon numbers, like the fluc- 
tuation of the net baryon number, take very different val- 
ues in the two phases of matter. To study dynamically 
the baryon number fluctuation in heavy ion collisions, 
we introduce a kinetic model that takes into account 
both production and annihilation of quark- ant iquark or 
baryon-antibaryon pairs. We first consider the case of 
baryon-antibaryon production from and annihilation to 
two mesons, i.e., •m\m,2 BB, in a hadronic matter of 
net baryon number s = Following the formalism 
of Ref. |l^] for the production of particles with conserv- 
ing charges, we have the following master equation for 
the multiplicity distribution of BB pairs: 

dP C 

^ = y(N mi )(N m2 ) (Pn-l-Pn) 
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In the above, P n {r) denotes the probability of finding n 
pairs of BB at time r; G = (ctgv) and L = {ulv) are 
the momentum-averaged cross sections for baryon pro- 
duction and annihilation, respectively; Nk represents the 
total number of particle species k; and V is the proper 
volume of the system. 

The equilibrium solution to Eq. (Q) is 



P„ 



"< oq 7 (2^)(^) 2 ' 
where I is the modified Bessel function, and 

c _ G(N mi )(N m2 ) 
L 

Using the generating function at equilibrium, 
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with g(l) = "^2 Pn = 1 due to normalization of the mul- 
tiplicity probability distribution, it is straightforward to 
obtain all moments of the equilibrium multiplicity distri- 
bution [fl3|| . For example, the mean baryon number per 
event is given by 



(.B)eq = b (n) = b 



Mil 

dx i 



bo^Y~ ~ b oy /e, (5) 
Jo 



In the above, bo is the fundamental unit of baryon num- 
ber in the matter, and the argument of Bessel functions 
J„'s is 2y/e. In obtaining the last expression in Eq. (||), 
we have kept only the leading term in y/e, correspond- 
ing to the grand canonical limit, y/e 3> 1, as baryons 
and antibaryons are abundantly produced in heavy ion 
collisions at RHIC 

The squared baryon number fluctuation is given by 
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It is seen that the mean number of baryons is propor- 
tional to bo while their squared fluctuation is proportional 
to Oq. The squared baryon number fluctuation per baryon 
is thus given by 



OJB,eq — bo 



bo 
2 ' 



(7) 



The mean number of antibaryons (B) and their fluctu- 
ation are the same as those of baryons as a result 
of zero net baryon number. Since bo is 1/3 in quark- 
gluon plasma and 1 in hadronic matter, wb and wg are 
smaller in the quark-gluon plasma than in an equilibrated 
hadronic matter by a factor of 3. 

We note that the equilibrium multiplicity distribu- 
tion in Eq. (|^) is not Poisson, as pointed out earlier in 
Rcf . jl6| . The non-Poisson distribution results from the 
quadratic dependence on the multiplicity n in the loss 
term of the master equation of Eq. (|l|) due to baryon 
number conservation. A Poisson distribution is obtained 
if the dependence on the multiplicity n is linear, which 
corresponds to production of particles that do not carry 
conserved charges. The master equation of Eq. ([l]) also 
gives a Poisson distribution at early times (r — » 0) when 
the loss term can be neglected. This corresponds to ei- 
ther production of particles with conserved charges dur- 
ing the early stage of heavy ion collisions or particle pro- 
duction without chemical equilibration as in e + e~ colli- 
sions. Since a Poisson multiplicity distribution gives 

D B ° issoa = b (B), (8) 

the squared baryon number fluctuation per baryon is 

Poisson i fn\ 



which is a factor of 2 larger than that in an equilibrated 
quark-gluon plasma or hadronic matter. A similar re- 
sult has been obtained previously by Gavin and Pruneau 
based on thermodynamic considerations jL7) . 



III. FINITE ACCEPTANCE AND NON-ZERO 
NET BARYON 

Because of experimental limitations, only protons and 
antiprotons in a certain rapidity and momentum range 
are usually measured. Moreover, the net baryon num- 
ber in heavy ion collisions is in general non-zero even at 
mid-rapidity due to the presence of projectile and target 
nucleons. To generalize the master equation to include 
these effects, we first consider the case of two species of 
antibaryons, e.g., antiproton and antineutron production 
from meson-meson interactions, i.e., myaii «-> Bp and 
m\mi <-> Bn in a hadronic matter of net baryon number 
s > [0. Defining P k ' 1 as the probability of finding 
k number of antibaryon species 1 and I number of an- 
tibaryon species 2 in an event (fc, I = 0, • • • , oc for s > 0), 
we then have the following generalized master equation: 



fc-iJ 
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dr ~ 

— — i [k(k +1 +s)P k - 1 - (k+l)(k + l +s+l)P k+1 ' 1 ] 
[l(k+l +s)P k ' 1 - {l + l)(k+l +s . 



(10) 



for 



In the above, 6j (j = 1,2) is defined as in Eq. 
the process involving antibaryon species j, and Gj and 
Lj denote the momentum-averaged cross sections for the 
gain and loss terms, respectively. We have neglected ef- 
fects due to quantum statistics, as they are not expected 
to be significant fill . 

The equilibrium solution to the above equation is given 
by the product of the equilibrium solution for the total 
BB pair number (k + l) and a binomial distribution, i.e., 



k.l 
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The distribution in Eq. ( |11[ ) can be understood intu- 
itively as resulting from first obtaining the distribution 
for the total antibaryon number (k + l), and then se- 
lecting k number of antibaryon species 1 and I number 
of antibaryon species 2 from the total antibaryons with 
probabilities of /i and fa, respectively. 

At equilibrium, the generating function 3(2/1,2/2) = 
Y J P k ' l y\y l 2 is g iven b Y : 
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9^y^)= {fiyi u2% J '( 2 vW/iyi+/altt)- (13) 

The above formulae can be equally applied to study 
the effects due to finite phase space coverage in experi- 
ments. In this case, antibaryon species 1 and 2 denote 
antibaryons in and out of the observed phase space, re- 
spectively; and /i and fi are their respective fractions. 

To generalize to N numbers of baryon and antibaryon 
species, we consider mimi <-> BiBj — 1,---,2N), 
where Bi and Bj with i.j — X,---,N denote baryon 
species i and antibaryon species j in the observed phase 
space, respectively, and Bi or Bj with i,j — N + 
1, • • • , 2N denote baryon species (i — N) and antibaryon 
species (j — N) out of the observed phase space, respec- 
tively. The joint generating function for both baryons 
and antibaryons is defined by 
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with the summation over all baryon and antibaryon num- 
bers, BiS and Bj's = 1, ■ ■ ■ ,2N). Assuming that 
particle momentum distributions are thermal, the gener- 
ating function at equilibrium is given by 



,V2n) 
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where 



with ejj defined for the process m,\mi «-> 
Eq. (|, and 



as in 



(17) 



From Eq. ([15j) , we obtain the mean numbers of baryons 
and antibaryons, 
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and their fluctuations, 
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In the above, we have kept only the first two leading 
terms in the expansions as s/^/e. <C 1 for ultra-relativistic 
heavy ion collisions. We note that Eqs. (|T^-p~9|) and ( pof - 
reduce to Eqs. (||) and (0) in the special case of 
full phase space coverage with only one baryon and an- 
tibaryon species (hi = fi = 1) and s = 0. 

Keeping only the leading terms in Eqs. ( pfl| ) and (|21 
leads to 



,eq 



1 



1 - 



(22) 



where hi and /j = 1, • • • , AT) represent, respectively, 
the fraction of total baryons that are observed baryon 
species i and the fraction of total antibaryons that are 
observed antibaryon species j. Eq.(p2[) thus shows that 
finite acceptance introduces a phase space correction fac- 
tor (1 — p/2), with p being the faction of the total baryon 
or antibaryon phase space that is observed, to the baryon 
or antibaryon number fluctuations in an equilibrated 
matter. We note that there is no such correction factor 
due to finite acceptance if the matter is far from equilib- 
rium when the BB multiplicity distribution is Poisson. 

FromEq. |[|), we can also derive the following squared 
net baryon number fluctuation for baryon species i and 
antibaryon species j: 



(16) (aN^) 2 ^((Bt-B,) 2 )-^-^ 



bofBi+Bj) 1 



(23) 



The squared net baryon fluctuation per baryon and an- 
tibaryon, (ANb) /(B + B), is thus a factor of 3 smaller 
in an equilibrated quark-gluon plasma than in an equi- 
librated hadronic matter, and is a useful observable as 
well H[ll|. If l ne P nase space acceptance is the same 
for baryons and antibaryons, i.e., hi = fj = p, then the 
phase space correction factor is (1 — p) as given in Ref. 



IV. CONCLUSIONS AND DISCUSSIONS 

We have shown that or cjg, the squared fluctua- 
tion of baryon or antibaryon numbers per baryon or an- 
tibaryon, reflects the fundamental units of baryon num- 
ber in the matter where they are produced. Since their 
expected values in the quark-gluon plasma are a factor 



3 



of 3 smaller than those in an equilibrated hadronic mat- 
ter, they can be used as signatures for the quark-gluon 
plasma formed in the initial stage of relativistic heavy 
ion collisions. Using a master equation with exact baryon 
number conservation, we have found that at equilibrium 
the baryon multiplicity distribution is not Poisson, and 
its squared fluctuation is half of that expected from a 
Poisson distribution. We have also calculated the cor- 
rection factors due to finite acceptance in experiments. 
Because of the contamination from the projectile and tar- 
get nucleons, lo§ is preferred over u>g. 

Although the baryon and antibaryon fluctuations are 
smaller in the quark-gluon plasma formed in the initial 
stage of heavy ion collisions, their magnitudes will be 
affected during the hadronization and subsequent scat- 
terings of baryons and antibaryons in the hadronic mat- 
ter. In the simplest quark coalescence model for the 
hadronization of the quark-gluon plasma to the hadronic 
matter, where three quarks form a baryon and a quark- 
antiquark pair forms a meson, the baryon number fluc- 
tuation is increased as a result of the finite fraction of 
quarks that form baryons. The correction factor to the 
baryon fluctuation in this hadronization scenario can be 
similarly derived as for finite acceptance. On the other 
hand, if hadronization is through the formation of strings 
and their subsequent fragmentation |20[| , then quark- 
antiquark pairs can be produced and annihilated dur- 
ing the string fragmentation and such additional pro- 
cesses usually also increase the baryon and antibaryon 
fluctuations. Furthermore, anomalous fluctuations could 
arise if the phase transition is first order [^1 22 ^or the 
quark matter is at the QCD tricritical point |7|23|]- As 
to the effects of hadronic scatterings, elastic scatterings 
of baryons and antibaryons are not expected to affect 
baryon and antibaryon number fluctuations in the full 
phase space. However, baryon-antibaryon pair produc- 
tion and annihilation tend to increase the fluctuations 
towards the values expected from a hadronic matter p4[ . 
The determination of baryon and antibaryon fluctuations 
in heavy ion collisions could be further complicated by 
the existence of gluonic baryon junctions p5|-p7|, so that 
baryon numbers are not necessarily carried fractionally 
by quarks but instead by the baryon junctions. In order 
for the baryon number fluctuation to be a viable signa- 
ture for the quark-gluon plasma, quantitative studies of 
above effects are needed. 
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